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The first-order phase transition in the three-state Potts 
model with long-range interactions decaying as l/r 1+CT has 
been examined by numerical simulations using recently pro- 
posed Luijten-Blote algorithm. By applying scaling argu- 
ments to the interface free energy, the Binder's fourth-order 
cumulant, and the specific heat maximum, the change in the 
character of the transition through variation of parameter a 
was studied. 



I. INTRODUCTION 

The critical behavior of models with long-range (LR) 
interactions has been considerably less explored than that 
of the short-range (SR) ones, especially for discrete mod- 
els where the non-locality of interactions makes most of 
the standard methods ineffective. These models may, 
however, exhibit rather complicated critical behavior al- 
ready in one dimension Recent attention has 
been driven to such models in the context of studying 
phenomena related to LR interactions but also in 
view of possible equivalence with SR models ||||] ■ 

One of the interesting and still not clarified aspects of 
LR models is the possible onset of the first-order tran- 
sition. The question naturally arises for the LR Potts 
model, which in the SR interaction case is known to un- 
dergo a first-order phase transition when the number of 
states q exceeds a certain limiting value which depends 
on dimensionality, q c (d) |p|-|lT|. 

We have recently pointed out that a similar be- 
havior indeed can be observed in the Id Potts model 
with interactions decaying with distance as l/r 1+<T and 
illustrated it on special cases q = 5, a = 0.2 and q = 
3, a = 0.8 characteristic of the two regimes. On the basis 
of these preliminary Monte Carlo (MC) results on small 
chains for q = 3 and q = 5, we have concluded on the ex- 
istence of a g-dependent threshold value a c (q) separating 
the first- and the second-order transition regime. 

In the present paper, we focus on a more detailed study 
of the Id three-state Potts model with < a < 1. 

As is well known, the first-order phase transitions are 
rather difficult to detect and study Most of stan- 

dard renormalization group (RG) approaches do not dis- 
tinguish them from the second-order transitions. 

On the other hand Monte Carlo simulations, which in 
combination with finite-size scaling (FSS) represent an 



efficient tool (m) to solve this problem, are limited when 
the first-order transition is weak, so that the correlation 
length, although finite, exceeds the considered system 
size. This problem is even more pronounced when deal- 
ing with LR interactions, where MC simulations become 
much more time consuming, no matter whether Metropo- 
lis or different cluster algorithms are used. 

Recently, Luijten and Blote [^5) have proposed a rather 
efficient algorithm which permits considering quite large 
sizes in spite of the long-range interactions. Our inten- 
tion is to use this approach here in order to perform 
a more systematic analysis of several quantities (inter- 
face free energy, Binder's fourth order cumulant, specific 
heat), which serve as criteria for distinguishing first- from 
second-order phase transition. 



II. MODEL AND METHOD 

The model we consider is defined by the Hamiltonian 
J 

Kj 



(1) 



where J > 0, Sj and Sj denote three-state Potts variables 
at sites i and j, respectively, S is the Kronecker symbol 
and summation is taken over all pairs of the system. 

This model has a phase transition at finite tempera- 
ture for < a < 1 No exact results exist for the 
related critical behavior. By analogy with the SR interac- 
tion case extended to arbitrary dimension, one expects to 
find a first-order phase transition for low values of a and a 
crossover to a second-order phase transition above a cer- 
tain threshold value a c . This has been confirmed by our 
preliminary MC results Jl2| in Id LR Potts model which 
lead to the conclusion, that for q > 2 the transition is of 
the first-order when a < <J c (q) and of second order above 
it. The threshold <J c (q) is expected to depend on q. For 
q = 3, the cr c was estimated [fl2| to lay above gmf = 0.5, 
the point separating the mean-field (MF) from the non- 
trivial critical regime in the Ising (q — 2) case ( where the 
transition, by symmetry reasons (see e.g. ref. J9| and ref- 
erences therein), remains of the second order fig ] in both 
regimes). Above the threshold a c , the model undergoes 
a second-order transition with non-classical critical ex- 
ponents depending on a. Several renormalization group 
approaches yield approximate values of these exponents, 
within continuous Ginzburg-Landau functional formal- 
ism |17j or in real space dJ|,0]. These latter approaches, 
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however, appear to be insensitive to detect a first-order 
transition. 

In the present approach we shall be using the MC al- 
gorithm, recently introduced by Luijten and Blote [ p"5| , 
designed for models with LR interactions. It is based on 
Wolf's cluster algorithm |Tj| and was applied to several 
MF related problems in the special (Ising) case q = 2 

IS 

The basic idea of the method is the use of cumulative 
bond probabilities in building the cluster of connected 
bonds, which drastically reduces the number of opera- 
tions required. The results coincide with those obtained 
with the simple Wolf algorithm, but the CPU time can 
be reduced by several orders of magnitude. 

The basic quantity in our calculations is the energy 
probability distribution defined by 



Pl(E) = 



1 



Zl(K) 



-KE 



(2) 



where K = l/T is the inverse temperature, J — kg, Z 
is the partition function, Nl{E) denotes the number of 
spin configurations corresponding to the energy E, and 
L is the system size. 

Several quantities used for the determination of the 
temperature driven first-order transition can be deduced 
from Pl(E). We concentrate here on three most impor- 
tant ones: the interface free energy, the Binder's fourth 
order cumulant poj , and specific heat. 

The interface free energy is obtained from the shape of 
the energy probabability distribution Pl(E). For tem- 
perature driven first-order transitions, at the transition 
temperature, Pl(E) has two maxima, corresponding to 
the coexisting ordered and disordered phases. The inter- 
face free energy is then defined by 



AF L = In 



Pl 



Max 



(3) 



where the finite-chain transition temperature K^ 1 has 
been defined by requiring that the two maxima are of 
equal height Plmux- Phmin denotes the minimum of 
Pl{E) between them. 

The scaling analysis of AFl can be used to identify 
the first-order phase transition, even in the case of a 
weak first-order transition, where the correlation length, 
though finite, is large and comparable to the system size 
pl| . When the transition is of the first-order, AFl in- 
creases with size. For systems with SR interactions the 
interface free energy has the dimension of surface and 
scales as AFl ~ L . In the present model with LR 
interactions it is expected to scale as a volume. 

The other two quantities can be derived from the 
higher energy momenta of Pl (E) , defined as 



(E n )L=J2 Enp L(E)- 



(4) 



The specific heat Cl (defined per spin in the remaining 
text) is related to the second moment and defined for the 
system size L by 



Cl = j?( (E 2 )l (E) 



(5) 



According to the FSS theory, in the case of a second-order 
phase transition its maximum scales as 



LMax ~ ) 



(6) 



while for the first-order one it simply scales as volume 
which in the present case means linear dependence on L 



ClMq 



L. 



(7) 



The calculation of the fourth moment gives Binder's 
fourth order cumulant [EOfl defined as 



V, 



(4) 



l-C/f/3 



where 



< E 2 >l 



(8) 



(9) 



In the present study, for practical reasons it is easier to 
deal with U L 4 \ In the thermodynamic limit U L ^ tends 
to one, when K ^ K c . At K — K c it still tends to one 
if the transition is of the second order, while it tends 
to a different constant in the case of a first-order transi- 
tion. Together with AF, appears as one of the most 
sensitive criteria for determination of first-order phase 
transitions lPQ£2l. 



A. The special case a = — 1 

Before proceeding to the presentation of our numerical 
results, let us summarize the only analytical results avail- 
able for the considered model. They can be obtained in 
the limit a — — 1 of model (1), where all the interactions 
are of equal strength and the coupling constant is rede- 
fined so that K — » K/L. This is the mean-field limit, 
which has been extensively studied in literature ( see e.g. 
Q ) and solved in the limit L — > oo by saddle point 
method || . 

The energy E and the entropy S of this model to the 
leading order in L writes 



E 
L 



1 



9-1 



2 ^ 

)»=() 



T 



S _ ^— \ L m j / L m 
L ~~ ~ ^ ~L n V ~ 

m=0 v 



(10) 



(11) 
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where L m 's, the numbers of particles in the m-th Potts 
state, satisfy the condition Lo + ... + = L. The 

transition temperature, the order parameter jump and 
the latent heat are known for the above model. 

To the best of our knowledge, the results for the char- 
acteristic quantities U$ ax and Cuax were not cited in 
literature. In order to obtain them, we have used the 
Lee and Kosterlitz [ p2[ prescription for the energy prob- 
ability distribution Pl (E) in the vicinity of the first-order 
transition temperature E~ x for large L, 



Pl{E) 



f (Lt) 5{E - E ) + 5(E - E do ) 
f(Lt) + 1 



(12) 



where t = (K — K c )/K c , and E Q and Edo are the ener- 
gies of ordered and disordered phases, respectively. The 
weight function f(Lt) which contains the entire tempera- 
ture and L dependence of Pl(E) is not given explicitly. 
Only the limiting behavior of f(Lt), f(x —* — oo) — > 
and f(x — > oo) — > +oo, is required, in order to have only 
one ^-function in thermodynamic limit for temperatures 
above or below the transition temperature, while for fi- 
nite L the two S- functions coexist. Maximizing U L and 
Cl with respect to temperature fixes the value of the 
function f(Lt) and leads to 



lim U 



(4) 

LMax 



(<Z-2) 4 



-i 2 



and 



lim 



Max 



2 (g-l) ( g 2 -3g + 3) 



■Mq-l) 



2q 



(13) 



(14) 



For finite L, the exact calculations may be performed 
numerically for very large sizes (L = 100 000 is easily 
attained). They match well with the analytical results 
( |l3| ) and ( |l4| ) and justify the conjecture Jl^ ) from which 
they were derived. 



III. NUMERICAL RESULTS 

By using the Luijten-Blote cluster algorithm, we have 
been able to reach sizes up to 5 000 with reasonable 
computing time. The systematic simulations were per- 
formed for sizes L ranging from 200 to 3 000, with peri- 
odic boundary conditions. The parameter a was taken 
with the increment of 0.1 in the interval < <r < 1, 
where the nontrivial transition is to be expected. For 
each set of parameters, 10 6 clusters were generated. Let 
us also point out that the characteristic temperatures re- 
lated to the three quantities considered are different un- 
til the thermodynamic limit is reached. Consequently, 
an independent numerical effort was needed for local- 
izing the characteristic temperatures for each of those 



quantities for every L and a considered. By combin- 
ing the Ferrenberg-Swendsen histogram method j24j] and 
direct calculations, these temperatures have been calcu- 
lated with numerical precision up to the fourth decimal 
digit in Kl , which implies a numerical error of approxi- 
mately 0.1%. We did not go beyond this precision in Kl 
since the aspect of the calculated distribution Pl{E) is 
much rougher than the one obtained with the Metropo- 
lis algorithm, when using a comparable number of steps. 
Resulting numerical error in calculated quantities them- 
selves varies from 1-5 % for Pl(E) and AFl to 1-2% for 

7 (4) 



A. Interface free energy 
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FIG. 1. Interface free energy divided by size versus inverse 
size. Data for a — 0.1 to 0.6 are included. Lines are obtained 
by linear regression. 



In Fig. 1 are presented the results for the interface free 
energy divided by size versus 1/L. Only sizes L > 400 
are included. Points of different shapes correspond to 
data with different a. It ranges from 0.1 to 0.6, where, 
within considered sizes, the two maxima of Pl(E) could 
be discerned beyond the error limits. For a = 0.6 this 
occurs only when L > 2 000. 

The lines represent the fit of AFl / L to linear form and 
illustrate the leading correction to the expected scaling 
form AF ~ L, and show good agreement with it. 

The weakening of the transition is manifested by the 
fact that AFl /L becomes smaller with increasing a. This 
dependence is almost linear unless a gets close to the 
assumed onset of the second order regime. 

The other two quantities, ujfi and Clmo,x, derived 



from momenta (Eqs. (|^) and (gj), have been calculated 
in the entire region < a < 1. 
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B. Binder's fourth order cumulant 



In order to find out whether in the thermodynamic 
limit U^Max wm tend to one or to a different value, we 
analyze the maxima of U^ ] {K) as functions of L. The 



1.15 




behavior is monotonous, the fit to the power-law form 



U 



(4) - c/ e (4) 



LMax 



const. ■ L x was made. The extrapolated 

results for ui^ are presented on Fig. 3. One can observe 
that the transition between the two regimes occurs in a 
continuous and smooth way with the change of a. For 
a < 0.6 we obtain ^ 1 outside the estimated error 
bars (given by the size of the points). Also, for a > 0.7 

(4) 

we conclude that, within the error bars, Ue has reached 
unity. The additional point 0.65 can be attributed to 
both regimes. 

These results are consistent with the above ones for 
AF. 
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FIG. 2. Fourth-order cumulant as a function of in- 
verse size. Points connected with dotted lines correspond 
to cr = 0.1,0.2,0.3,0.4,0.5,0.6,0.65,0.7,0.8,0.9, from top to 
bottom. The solid line denotes the results for a = — 1 limit. 



results for U 



(4) 

LMax 



are summarized in Fig. 2, where the 



a = — 1 results are also traced for comparison. For low 
values of a, U L ^ Max is non-monotonic in 1/L, similar to 
the a = — 1 case, so that in the L — > oo limit it is clearly 
different from unity. For higher values of cr, where the 
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C. Specific heat 

The data for the specific heat maxima C^Max are sum- 
marized in Fig. 4 in the form of C'lmux/L versus 1/L. 
According to Eqs. M) and (R) and by taking into ac- 
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FIG. 4. Maxima of the specific heat divided by 

size versus inverse size. Solid line denotes the case 
a = —1. Points connected with dotted lines correspond to 
a = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.65, 0.7, 0.8, 0.9, from top to bot- 
tom. 
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FIG. 3. Extrapolated values of the fourth-order cumulant 
as a function of a. The constant line at unity is drawn to 
guide the eye. 



count that a/v is much smaller than unity in the region 
of interest, ChMax/L should tend to a constant or to zero 
value, depending on whether the transition is of the first- 
or second-order, respectively. 

Similar as for Uj/, the curves are non-monotonic and 
tend to a nonzero value in the a = — 1 case and for 
low values of a. For higher values of a the convergence 
to limit L — > oo was examined by the fit to the form 
ChMax/L — Co + const. ■ L~ x . The constant cq is nonzero 
in the first-order regime, while in the second-order one 
co = and x = 1 — a/v. This means also that the 
log-log plot of the curves of Fig. 4 should have a linear 
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shape in the latter case, while it will deviate from the 
straight line in the former case. In Fig. 5 log-log plots of 
Clmo,x/L versus 1/L are given for a > 0.5. The curves 
for 0.5 and 0.6 clearly leave the straight line. Linear 
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FIG. 5. Log-log plot of the maxima of the specific heat 
versus inverse size. The points connected by solid lines rep- 
resent the results for a = 0.5,0.6,0.65,0.7,0.8 and 0.9, from 
top to bottom. Dotted lines denote the fit to the straight line 
by linear regression. 

fit is very good for a > 0.65, indicating the second- 
order transition. The corresponding critical exponent 
ajv following from the fit for a = 0.65,0.7,0.8,.0.9 is 
0.36,0.33,0.26,0.19, respectively, which is in poor agree- 
ment with values expected from the finite-range scaling 
(FRS) H calculations, which for the same values of a 
would give ajv = 0.32,0.27,0.15,-0.02. One should 
mention, however, that it is generally a difficult task to 
extract critical exponents with good precision from the 
spacific heat maxima p6[ . 

D. Critical temperature 

The three different characteristic temperatures arise in 
calculation of AFl and maxima of U^ ] {K) and C L (K), 
which all should tend to the critical temperature in the 
thermodynamic limit. The L — » oo extrapolations of 
these quantities made by assuming power-law corrections 
in (1/L) are presented in Table I, compared to the ex- 
isting earlier results obtained by FRS [^) and RG 0. 
The extrapolation errors are estimated to be in the limits 
(±0.01) which is also comparable to difference between 
the extrapolations obtained from different characteristic 
temperatures. The agreement with the FRS results is 
quite good, discrepancy does not exceed 5 %, except 
for a = 0.1, where the convergence of FRS is weakest, 
and the extrapolation errors in both methods are largest. 



(For a = 0.1 all the claculations were thus performend 
up to the sizes L = 5000.) The RG results appear to 
be systematically larger then ours, and the discrepancy 
with them varies from 7 to 33 %. Similar discrepancy is 
obtained for the fit to the functional form K c ~ a in the 
limit a — ► 0, conjectured in ref. ||. 

TABLE I. Inverse critical temperatures (K e (C), K e (U {4) ), 
K e (AF)), obtained by extrapolation of Kl compared to FRS 
extrapolated values (K(FRS)) || and RG results (K(CM)) 



K e {C) K e (U w ) K e (AF) K(FRS) K e (CM) 



0.1 


0.16 


0.15 


0.16 


0.14 


0.15 


0.2 


0.28 


0.28 


0.27 


0.270 




0.3 


0.38 


0.38 


0.37 


0.386 


0.43 


0.4 


0.49 


0.49 


0.48 


0.494 




0.5 


0.58 


0.60 


0.59 


0.601 


0.71 


0.6 


0.71 


0.71 


0.71 


0.714 




0.65 


0.78 


0.77 




0.774 




0.7 


0.84 


0.84 




0.837 


1.05 


0.8 


0.98 


0.99 




0.977 




0.9 


1.14 


1.13 




1.144 


1.64 



IV. CONCLUSION 

We have applied MC simulations in combination with 
FSS in order to examine the onset of the first-order phase 
transition in the Id three-state Potts model with long- 
range interactions. The Luijtcn-Blote advanced algo- 
rithm for long-range interaction systems allowed us to 
treat successfully considerably large sizes (up to 3 000) in 
a reasonable amount of time, in spite of the long range of 
interactions. The systematic analysis of three quantities: 
the interface free energy, Binder's fourth order cumulant, 
and specific heat, confirms the existence of two regimes 
in the interval < a < 1 . All three considered quantities 
give the first-order transition for a < 0.6. The transition 
becomes gradually weaker with increasing a and changes 
smoothly to a second-order transition by continuous vari- 
ation of a. 

For the sizes examined here (up to 3 000) and within 
the estimated error bars, we obtain a second-order tran- 
sition for a > 0.7. 

The present approach has, however, limitations to pre- 
cise determination of the threshold. Since the transition 
close to the threshold becomes arbitrarily weak, the finite 
correlation length, characteristic of the first-order tran- 
sition, will always be much larger than the size of the 



considered system, close enough to a c . Thus, the present 
result for cr c should be understood only as a lower limit 
for the possible onset of the second-order phase transi- 
tion. 

For this purpose rather complementary studies should 
be done, such as the one on the dependence of the finite 
correlation length on a. 

It is interesting to notice at the end that, according 
to present results, o~ c {q — 3) falls in the interval between 
0.6 and 0.7. If the correspondence between SR and LR 
models pi would be extended outside the MF regime to 
the present problem and to the Potts model, it would 
lead to the conjecture o c {q = 3) = 2/d c (q = 3)sr which 
would give a c (q — 3) close to and slightly larger than 
0.66. However, this line of argument would also imply 
that <7 C = 1 already for q = 4. 
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